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We experimentally examine a laser-induced underwater shock wave with a special atten¬ 
tion to pressure impulse, the time integral of pressure evolution. Plasma formation, shock- 
wave expansion, and pressure in water are observed simultaneously using a combined 
measurement system that obtains high-resolution nanosecond-order image sequences. 
These detailed measurements reveal a non-spherically-symmetric distribution of pressure 
peak. In contrast, remarkably, pressure impulse is found to distribute symmetrically for 
a wide range of experimental parameters even when the shock waves are emitted from 
an elongated region. The structure is determined to be a collection of multiple spherical 
shock waves originated from point-like plasmas in the elongated region. 
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1. Introduction 

Underwater shock waves induced by illumination with a nanosecond laser pulse are 


utilized in various applications including low-invasive medical treatments ([Razvi et al 


1996| Sofer et a/.||2002t Lam et a/.||2002| |Sankin et a/.||2005| Klaseboer et a/.||2007t |Lee & 

Doukas 1999| Kodama et a/.||2000). The laser-induced shock wave ca 

n trigger a sudden 
ation of high-speed 

motion of a liquid on a free surface, which, for instance, results in gener 

microjets applicable for needle-free injection devices (Menezes et al. 

2009 

Thoroddsen 


2013 Marston fc Thoroddsen||2Q15). 


For the sudden motion of the liquid, one of the most important quantities is pressure 
impulse ( Batchelor|1967 Cooker fc Peregrine|1995 Antkowiak et a/.|2QQ7 ). Its definition 
is given as: 


-I 


P= pdt, 


( 1 . 1 ) 


where p is pressure of the liquid and t is the elapsed time. Peters et al. (2013) numerically 


reproduced the high-speed microjet reported by Tagawa et al. (2012) and confirmed that 


the pressure impulse is the key quantity for the motion of the jet. Thus detailed investi¬ 
gation for the pressure impulse of the laser-induced shock wave is of great importance. 

The shock wave has been often modeled as a spherieal shoek, which assumes a spherically- 
symmetric pressure distribution and spherical shape of the shock. However, some re¬ 
searchers have pointed out that the spherical-shock model is not applicable in cer¬ 


tain cases (Buzukov et al. 1969 Sankin et al. 2008 Noack & Vogel 1998 Vogel et al 
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Figure 1. Snapshots of (i) a collection of plasmas in a conical region with 10 x objective, 6.9 
mJ and (ii) a collection of shock waves originated from plasmas. Time t denotes the elapsed 
time from the start of illumination with the laser pulse. 


1996a Lauterborn fc VogelpQld ). Buzukov et al. (1969) reported that a non-spherically- 


symmetric bubble is observed with a series of compression waves. Sankin et al. (2008) 
measured pressure peaks for a shock at various positions and determined that the peak 
pressure at a point in the direction perpendicular to the laser beam is more than twice 


as high as that in the direction of the laser. Vogel et al. (1996a); Noack & Vogel (1998) 


reported that the shape of a shock wave is not spherical due to conical plasma forma¬ 
tion. The movements of the breakdown front during plasma formation in water had been 
intensively studied utilizing streak photographs (e.g. Docchio et al. |( |19886|a| | ).In our ex¬ 
periments we also observe a non-spherically-symmetric plasma, bubbles, and shock waves 
as shown in Figure Figure [^i) shows a collection of point-like plasmas in a conical re¬ 
gion (like “grapes of plasmas”). Figure [^ii) shows another snapshot of both bubbles and 
shock waves. It confirms that for small numerical apertures (NAs) of an optical system 
the shock waves are not from a single conical plasma but from a collection of plasmas in 
a conical region. Note that Vogel et al. (19966) did show single conical plasmas. Despite 
these considerations, a common model for the pressure impulse of laser-induced shock 
waves has not been developed. 

In this study, we report on experimental observations of a laser-induced shock wave 
with a special attention to pressure impulse. We also propose a new model of the shock 
wave to rationalize the observations. Such an observation is, however, challenging because 
each phenomenon involved in generating the shock occurs within a short time; illumina¬ 
tion with a laser pulse first triggers the emergence of plasma in water, which leads to 
rapid expansion of a bubble and emission of a shock wave ( Noack fc Vogel||1999 


Lauter- 


born et ^|2QQ1 ). The time scale for plasma growth is in the order of nanoseconds and 


the shock velocity in water is approximately 1,500 m/s. In this study plasma growth, the 
expansion process of the shock, and pressure in water have been simultaneously measured 
using a combined measurement system, in which ultra-high-speed recording systems and 
pressure sensors are installed. 


2. Experimental setup and method 

Figure shows the combined measurement system. An underwater shock wave is in¬ 
duced by a 532 nm, 6 ns laser pulse (NdiYAG laser Nano S PIV, Litron Lasers) focused 
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through an objective lens to a point inside a water-filled glass container (100x100x450 
mm). The initial laser beam diameter is 4 mm. Water is distilled by a water-purification 
system (Milli-Q Integral, Merck) at room temperature (15^20 °C) and gas-saturated. 
Its electrical conductivity is 13 MQ- cm. The two experimental parameters are the mag¬ 
nification of objective lens (5x [NA 0.1], 10x [NA 0.25], 20x [NA 0.25], MPLN series, 
Olympus) and the laser energy (2.6 mJ, 6.9 mJ, 12.3 mJ). The parfocalizing distance of 
the objective (PFD) is 45 mm for all the microscope objectives while working distances 
are 20 mm, 10.6 mm, and 25 mm for 5x, 10x, and 20x objectives, respectively. Focus¬ 
ing angles of each microscope objective are 1 degree, 4 degrees, and 6 degrees for 5x, 
10 X, and 20 x objectives, respectively. The diffraction-limited focused-beam diameter d 
is calculated by following equation (Vogel et al. ( 2005[ )), 

A 


d= 1.22 


NA' 


( 2 . 1 ) 


where A (= 532 nm) is the wavelength of a laser, NA (5x [0.10], 10x [0.25], 20x [0.25]) 
is the numerical aperture of microscope objective. The calculated d is 6.5 /im, 2.6 /im, 
and 2.6 jim with 5x, 10x, and 20x, respectively. 

The combined measurement system consists of two hydrophones (Muller Platte-Gauge, 
Muller) and two ultra-high-speed cameras. One of the hydrophones is placed 5 mm away 
from the focal point of the laser in the direction of the laser beam {0 = 0° direction). The 
other hydrophone is at the same distance but in the direction perpendicular to the laser 
beam {0 = 90° direction). The impulse response time (rise time of an impulsive signal) of 
the hydrophones (the piezoelectric PVDF type hydrophone) utilized in this study is 35- 
45 ns. The hydrophones are connected to an oscilloscope (ViewGo II DS-5554A, Iwatsu) 
for recording hydrophones’ signals. The sampling frequency of the oscilloscope is 2 GHz, 
temporal resolution of 0.5 ns. The oscilloscope digitalizes pressure value every 0.03 MPa. 
One of the cameras is an ultra-high-speed camera (Imacon 200, DRS Hadland) with up to 
200x10^ fps (5 ns time interval) and a 1,200x980 pixel array to record plasma formation 
and shock waves in the near field. The other camera is another ultra-high-speed video 
camera (Kirana, Specialized Imaging) with up to 5x10^ fps and a 924x768 pixel array 
for imaging shadowgraph of shock-wave propagation. This camera is synchronized with 
a laser stroboscope that operates with a pulse width of 20 ns as a back illumination 
source (SI-LUX 640, Specialized Imaging), the repetition rate of which is also up to 
5x10^ Hz. A digital delay generator (Model 575, BNG) is used to trigger the laser, the 
hydrophones, the cameras, and the stroboscope. Each measurement was repeated more 
than three times under the same experimental conditions. 


3. Results and discussion 

3.1. Observations and pressure measurement 

Figure]^ (in which t denotes the elapsed time from the start of illumination with the laser 
pulse) shows the measurement results obtained with the 10 x objective lens. Figure [^i) 
shows a snapshot of the plasma luminescence in an elongated area, the major axis of 
which is in the direction of the laser beam. The image sequence of the plasma confirms 
that all parts of the plasma emit strong lights within ±5 ns. A laser-induced bubble 
then emerges where the plasma was formed (Fig.|^ii)) and its shape is also elongated in 
the direction of the laser beam. At t = 0.4 /is, non-single spherical shocks are observed. 
In contrast, at t = 2.4 /is, the shock could be regarded as a single spherical shock (see 
Fig. [^iii)). However, enlarged images for 6> = 0° and 0 = 90° (Fig, [^ih) 6> = 0° and 0 
= 90°) display a clear difference. Two shock waves for ^ = 0° (Fig. I^ih) 0 = 0°), which 
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Ultra-high-speed video camera (200 Mfps) 


Figure 2. Measurement system consisting of two ultra-high speed cameras and two pressure 
sensors. An ultra-high speed video camera records laser-induced shock waves and bubbles at up 
to 5x10® fps with a synchronized laser stroboscope. Plasma luminescence is captured by another 
ultra-high speed video camera at up to 200x10® fps. Temporal pressure evolution is measured 
by two hydrophones. One hydrophone is arranged in the direction of the laser beam (0 = 0°) 
at a stand-off distance of ca. 5 mm from the laser focal point, while the other hydrophone is at 
right angles to the hydrophone (0 = 90°) at the same stand-off distance. 


is different from the single shock wave for 0 = 90° (Fig. [^ih) 0 = 90°). Figure [^iv) 
shows the temporal evolution of pressure measured with the two hydrophones placed at 
different positions. There are two peaks for 0 = 0°, while there is a single large peak for 
0 = 90°, which is approximately 1.3 times higher than that for ^ = 0°. Note that this 
dependence of the peak pressure on the angle 0 is the same as that reported by |Sankm| 
et al. (2008). For the 5x objective lens, this trend is much more pronounced: four plasma 
groups and four shock waves separated from each other are evident (see Fig. [^. 


3.2. Pressure impulse 

Here, we compute the pressure impulse for 0 = 0°, Pq, and that for 0 = 90°, Pgg. 
We calculate pressure impulse from t = 2.5 to 4.5 //s, during which time-lags of plasma 
formation (<10 ns) is totally covered, i.e. the time-lags do not affect the results of pressure 
impulse. Integrations for the pressure with respect to the elapsed time indicate pressure 
impulses. Both Pq and Pgo for the shock obtained with 10 x objective lens are shown 
in Fig. [^iv) in the unit of Pa-s. Pg is in reasonable agreement with Pgo- Furthermore, 
pressure impulse and peak pressure for 0 = 0° and 90° were examined for all the other 
experimental conditions. While peak pressure for ^ = 0° and 90° differ significantly 
as shown in Fig|^i), the Pg was in agreement with the corresponding Pgg within the 
experimental uncertainty for a wide range of experimental parameters (Figj^ii)). Note 
that the order of pressure impulse in this study is the same order of the practical use for 
drug delivery for cytoplasmic molecules ( Kodama et a/.|20QQ ) and generation of micro jets 
( Tagawa et a/. [20 12 2013). 
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t=5ns f=0.4|js ^ = 2.4|JS 



Time, |js 


Figure 3. Measurement results for a laser-induced underwater shock wave obtained with a lOx 
objective lens, (i) Plasma luminescence at t = 5 ns after the laser is fired. The image is captured 
by an ultra-high-speed video camera at 200 Mfps. (ii) Shock waves and bubbles at t = 0.4 //s 
imaged with an ultra-high speed video camera at 5 Mfps. (iii) Shock waves at t = 2.4 /xs measured 
with an ultra-high-speed video camera at 5 Mfps. Enlarged images for the areas of ^ = 0° and 
0 = 90° are also presented, (iv) Shock pressure measured by the two hydrophones arranged at 
^ = 0° (red line) and 0 = 90° (blue line) with respect to the laser direction. Integrations for the 
pressure with respect to the elapsed time indicate pressure impulses. 


3.3. Strueture of multiple shoek waves 

Based on the aforementioned results, we here propose a model for the structure of the 
laser-induced shock wave: The shock has a multiple structure that consists of multiple 
spherical shock waves as depicted schematically in Fig. We assume that each spherical 
shock wave originates from the corresponding plasma formation. In addition, since the 
shock wave behaves acoustically with low pressure (< 100 MPa ( Vogel et a/.|1996^ ), we 
could apply the superposition prineiple to analyze the pressure impulse in the far-held. 
The net pressure that is produced by two or more shock waves reaching the same point 
is the sum of the pressure induced by the individual shock waves. A phenomenon that 
is analogous to this may be the surface wave observed after one or several stones are 
thrown into a quiescent pond (the so-called Huygens-Fresnel principle). Note that the 
Huygens-Fresnel principle does not apply for nonlinear shock wave propagation. 

For both the 10 x and 5x objective lens, this model rationalizes the observations of 
both pressure peaks and pressure impulse: several peaks for 6> = 0° and the single large 
peak for 0 = 90° while pressure impulse for 6> = 0° matches pressure impulse for 0 = 90°. 
Results for a wide range of experimental parameters (Fig[^ may indicate the universality 
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(i) (ii) (iii) 



t=5r\s ^=0.4|js t = 2A[}s 

(iv) 



Time, [}s 


Figure 4. Measurement results for a laser-induced underwater shock wave obtained with a 
5x objective lens. Captions for (i)-(iv) are the same as those in Fig.[^ 


of the multiple structure model for optical breakdowu at low or moderate NA. Note that 
this sceuario iucludes the well-kuowu spherical-shock model. It should be emphasized 
that, as observed with the 10x objective leus (Fig. [^, eveu if just a siugle plasma or 
a bubble is observed, the origiu of the elougated plasma aud bubble is expected to be 
multiple spots of plasma, which leads to the emergeuce of multiple spherical shocks, 
resultiug iu a uotable augular variatiou of shock pressure. This model for the multiple 
shock structure could possibly ratioualize the results reported iu previous research. For 
iustauce, Saukiu et al (2008) reported that a laser-iuduced shock wave emitted from 
au elougated plasma has au augular variatiou of pressure distributiou iu the far filed. 
Although the shape of the shock wave appears to be spherical, the elougated plasma 
may cause a multiple-structure of the shock, as observed iu the preseut experimeuts (see 
Fig. §, which would lead to a uou-spherically-symmetric pressure peak of the shock. 
Besides aforemeutioued pheuomeua iu the far field, au auisotropy of shock wave iu the 
uear field of cyliudrical plasmas had beeu reported (e.g. Schoeffmauu et al ( 1988[ )), which 
we discuss iu detail iu Sec. ' ' 


3.4. Laser-induced plasma and bubbles 

Here, atteutiou is giveu to plasma formatiou, which is the origiu of the multiple shock 
structure. Figure [^i) shows the plasma lumiuesceuce uuder all the experimeutal coudi- 
tious. Figure [^ii) shows the leugth of the plasma regiou iu the directiou to the laser beam 
as a fuuctiou of laser euergy for 5x, 10x, aud 20x objectives. The leugth iucreases with 
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5x objective lens 10x objective lens 20x objective lens 

Figure 5. (i) Peak pressure of a laser-induced underwater shock wave measured at ^ = 0° and 0 
= 90° for all the experimental conditions. Each presented value is a mean for three measurements 
and its error bar shows the standard deviation, (ii) Pressure impulse of a laser-induced shock 
wave measured at ^ = 0° and 0 = 90° for all the experimental conditions. 


the laser energy for all objectives. With lower magnification, the length of the plasma 
region becomes longer. The relation between an elongated plasma length (region) and 
focusing angle or laser pulse energy had been intensively investigated experimentally by 


Vogel et al. (19966). They had found that there is a strong dependence of the elongated 


plasma length on the focusing angle and laser energy: Plasma is created in the cone an¬ 
gle of the laser beam proximal to the laser and becomes longer with larger laser energy. 
Our results are consistent with theirs. Thus the multiple plasma formation is dependent 
on the spherical aberrations of the focusing optics, liquid impurities and the input laser 
energy ( Vogel et a/.||1996a|6 Nahen fc Vogel|p~996 ). 

There are two mechanisms which can lead to plasma formation: direct ionization of the 
medium by multiphoton absorption or avalanche ionization via inverse bremsstrahlung 
absorption ( [Vogel et~aT. (19966)). The mechanism leading to optical breakdown (plasma 
formation) in this study is the avalanche ionization by heating of impurities since the 
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Figure 6. Schematic of a multiple structure model for a laser-induced underwater shock wave. 
Multiple plasmas emit multiple bubbles and spherical shock waves. The shock pressure at a 
certain point is the sum of spherical shock pressures that reach the same point. 




Laser energy [mJ] 


Figure 7. (i) Plasma formation under all experimental conditions (Magnihcation of objective 
lens: 5x, 10x, 20x. Laser energy: 2.6, 6.9, and 12.3 mJ). The plasma shape is the most elongated 
with an input energy of 12.3 mJ and the 5x objective lens, whereas it is rather spherical with 
an input energy of 2.6 mJ and the 20x objective lens, (ii) Length of plasma as a function of 
the laser energy. The circle plot and error bar show respectively the mean and the standard 
deviation in 5 trials. Red, blue and green color bar represent the microscope objective of 5x, 
10 X and 20 x, respectively. 
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(i) 5x (ii) 10x 



Figure 8. Snapshots of laser-induced bubbles and shock waves with 5x, lOx, 6.9 mJ, t = 50 
ns. These results are obtained under the same condition. 


position of plasma formation varies within a certain region. Note that impurities may 
provide centers for both linear and nonlinear absorption. If the multiphoton process 
in nonlinear absorption has a lower order than in pure water, impurities can trigger 
breakdown. Figure shows snapshots of bubbles (and shock waves) under the same 
condition (the same magnification, laser energy, etc.). The number/position of plasma 
vary at every laser shot. If the plasma formation mechanism were purely direct ionization 
of the medium by multiphoton absorption, identical plasma shape should be observed 
for every laser shot in the same focusing optics with the same laser energy. 

Note that the multiple plasma generation will not be avoided by the “perfect” focusing 
without any spherical aberrations since the plasma occurs in the region where the local 
energy exceeds the breakdown thresholds and thus the plasmas do not always occur at 
the perfectly focused point (Vogel et al. 1999a 2005). Therefore a focusing angle is a 
crucial parameter since it strongly affects local energy density. In our experiment, there 
exists a strong relation between the length of the plasma region and focusing angle. 

In the case of linear sound propagation, for an extremely elongated cylindrical (or 
conical) source, anisotropic emission is expected with more total energy and total impulse 
in the directions perpendicular to the cylinder than along its axis. However both bubbles 
and shock waves in the near field (e.g. Fig.[^ show that the source is not a single plasma 
but a collection of point-like plasmas in a conical region. Therefore, even in the case of 
linear sound propagation, it is not quite obvious to assume anisotropy of total energy 
and impulse. 

The number of bubbles depends on the laser energy and the objective magnification as 
shown in Fig.[^ The number of bubbles increases with increasing laser energy and lower¬ 
ing the objective magnification. This trend is consistent with the length of plasma region 
as shown in Fig. We measure distances between bubble centers and count number of 
times that a certain distance emerges for all experimental conditions. Figure [T^a), (b) 
and (c) represent the histogram of the distances of each bubble centers for the microscope 
objective of 5x, 10x and 20x, respectively. The number of times is the sum of three tri¬ 
als. The average distance for each condition depend on the focusing angle (microscope 
objective) and the laser energy. Bubble centers are not always equidistant. 
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2.6 mJ 6.9 mJ 12.2 mJ 
Laser energy [mJ] 


Figure 9. Number of bubbles of as a function of the laser energy. Each color bar represents a 
result of one pulsed laser. Red, blue and green color bar represent the microscope objective of 
5x, 10X and 20x, respectively. 



Figure 10. The histogram of the distances between bubble centers. Vertical axis is number of 
times that a certain distance emerges for all experiments, (a), (b) and (c) show the histogram 
of the microscope objective of 5x, 10x and 20x, respectively. The color of bar represents the 
energy of a pulsed laser. The blue, green and red bar mean 2.6, 6.9 and 12.3 mJ. 


3.5. Shock wave in the near field 

In this section, we estimate the pressure decay in the near field to discuss energy dis¬ 
sipation. The estimation of energy dissipation for multiple plasmas and shock waves 
is, however, a big challenge since it is hard to disentangle the energy dissipation and 
nonlinear interaction in the near field. 

We now use a high-speed camera for following shock wave emission from the source in 
the near field. We take snapshots of the shock wave, follow its position (Fig.pTj^i)), derive 
its velocity (Fig. [T^ii)), and estimate the pres sure (Fig. pTj ^iii)). To estimate pressure in 
a shock front Ps, we apply following equation (Vogel et al. ( |1996a| )). 


Ps = CiPqUs{\Q^^‘ - 1) +Poo, 


(3.1) 
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where po is the density of water before compression by the shock wave, cq is the normal 
sound velocity in water, ci = 5,190 m/s, C 2 = 25,305 m/s, and p^^ is the hydrostatic 
pressure. The results show that, in the near field, pressure of the shock in the direction 
perpendicular to the laser beam is higher than that in the direction of the laser beam as 


shown in Fig.llUiii). Note that Vogel et al (1996a); Noack & Vogel (1998) found similar 


pressure values close to the plasma rim as presented in the present study. The energy 
dissipation at the fronts of high-pressure shock waves is thus expected to be faster than 


at the low-pressure fronts since Vogel et al (19996) showed that the dissipation rate 
of acoustic energy is proportional to the pressure jump at the shock front. However the 
pressure decay related to energy dissipation is slightly slower at the fronts of high-pressure 
shock wave than that at the low-pressure fronts. Similar results had been reported by 


Schoeffmann et al. (1988); Vogel et al. (1996 a) that the pressure decay is not faster at the 


fronts of high-pressure shock wave than that at the low-pressure fronts in a certain case. 


Vogel et al. (1996 a) attributed it to the formation process of the shock wave: the pressure 


maximum is located behind the leading edge of the pressure transient. Our interpretation 
for this is that, in the near field along the direction perpendicular to the laser beam, a 
shock wave from one of the plasmas is followed by the other shock waves originated from 
the other plasmas that eventually overlap and add pressure of the shock with increasing 
the distance. This might lead to the slower pressure decay compared to the case of a 
single shock wave. 

It is noteworthy that high pressures in the near field combined with nonlinear propa¬ 
gation and a strong anisotropy translate into the linear isotropic pressure impulse in the 
far field. We discuss mechanisms of this interesting finding in this paragraph. For a cylin¬ 


drical source Schoeffmann et al. (1988) showed an anisotropic shock wave emission with 


most energy in the 90° direction due to the geometrical effect. The anisotropic pressure 
jump at the initial shock front suggests that the total amount of energy dissipation is 
higher in the 90° direction than in the 0° direction (cf. [Vogel et al. (19996)). Therefore 
one likely explanation for the transition is that the initial anisotropy is eroded with in¬ 
creasing propagation distance during nonlinear sound propagation in the transition from 
near field to far field, and an isotropic impulse distribution in the far field can evolve. 
Another possible explanation is that nonlinear propagation in the near field has limited 
effects on the isotropic impulse distribution in the far field under the present experi¬ 
mental conditions. In order to address aforementioned discussions, further experimental 
evidence of ultra-high speed recordings would be needed. 


4. Conclusion and Outlook 

In order to investigate a laser-induced shock wave with focusing on pressure impulse, we 
constructed a measurement system consisting of a combination of ultra-high-speed cam¬ 
eras and pressure sensors. Shock pressure was measured with two hydrophones arranged 
at ^ = 0° and 0 = 90° with respect to the laser direction, and plasma formation and shock 
wave expansion were simultaneously observed using two ultra-high-speed video cameras. 
The most important finding in this paper is that the distribution of pressure impulse 
of a shock wave is spherically symmetric (isotropic) for a wide range of experimental 
parameters even when the distribution of peak pressure is non-spherically-symmetric 
(anisotropic). We proposed a multiple structure model for the laser-induced shock wave: 
The laser-induced shock wave is a collection of spherical shock waves emitted from the 
bright spots inside the area of plasma luminescence. The multiple structure is dependent 
on the plasma shape generated by illumination with the laser pulse. 

To the best of the authors’ knowledge, the isotropy of the pressure impulse in the far 
































12 


140 


120 


V. Tagawa, S. Yamamoto, K. Hayasaka, and M. Kameda 

(iii) 


i 100 


0 

o 

S 80 

CO 


40 



10 ^ 


CO 

CL 


0 

□I 


— 0 = 0 ° 
— 0 = 90 ° 



102 L 
101 


102 103 

Distance [|jnn] 


Figure 11. (i) The distance between plasma and a shock front vs. the elapsed time, (ii) The 
velocity of the shock wave vs. the elapsed time, (iii) The pressure of the shock wave vs. the 
distance between plasma and a shock front. The inserted numbers indicate the local slope of the 
corresponding curve. 

field is reported for the first time, which is of great importance for various applications. 
For instance, the pressure impulse in this study is in the same order of the practical use 
for drug delivery for cytoplasmic molecules (Kodama et al (2000)). Other examples are 


low-invasive medical treatments such as drug injection and lithotripsy, for which the in¬ 
put energy of the laser is of the order of 10 mJ (Tagawa et al ( 2012[[^13| )) or even more 
(Menezes et al. (2009)). The laser energy in this study is in the same range: 2.6 - 12.3 mJ. 


The isotropic distribution of pressure impulse may provide high degrees of freedom for 
the design of needle-free injection devices using high-speed microjets. By changing the 
plasma shape with control parameters (magnification of the objective lens or the input 
laser energy), the anisotropy of the shock pressure could be controlled, which might be 
applicable to a variety of advanced techniques. 
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